logist_binom_predict <- posterior_predict(stan_logistic_binom, draws = 500)
str(logist_binom_predict)
## 'ppd' int [1:500, 1:19] 1231 1220 1196 1224 1245 1194 1261 1226 1224 1226 ...
## - attr(*, "dimnames")=List of 2
## ..$ : NULL
## ..$ : chr [1:19] "1" "2" "3" "4" ...
# convert sucesses to proportions
logist_binom_predict <- t(t(logist_binom_predict) / golf$tries)
str(logist_binom_predict)
## ppd [1:500, 1:19] 0.853 0.845 0.829 0.848 0.863 ...
## - attr(*, "dimnames")=List of 2
## ..$ : NULL
## ..$ : chr [1:19] "1" "2" "3" "4" ...